The Superconducting Quantum Point Contact 
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I. INTRODUCTION 

The work reviewed in this contribution was motivated 
by a sequence of analogies. First came the quantum point 
contact (QPC) for electrons 0,0, a constriction in a two- 
dimensional electron gas with a quantized conductance. 
Then followed the optical analogue: The predicted Q dis- 
cretization of the transmission cross-section of a diffusely 
illuminated aperture or slit was recently observed exper- 
imentally |^. One can speak of this optical analogue as 
a QPC for photons. Can one extend the analogue to- 
wards a QPC for Cooper pairs? The answer is Yes p|- A 
narrow and short, impurity- free constriction in a super- 
conductor has a zero-temperature critical current which 
is an integer multiple of eAo/h, with Ao the energy gap 
of the bulk superconductor. In this paper we present our 
theory for the discretization of the supercurrent, in more 
detail than we could in our original publication [5( . For a 
less detailed, but more tutorial presentation, see Ref. 0. 
To introduce the reader to the problem, we first summa- 
rize some basic facts on the Josephson effect in classical 
point contacts 0. 

The theory of the stationary Josephson effect for a 
classical ballistic point contact is due to Kulik and 
Omel'yanchuk The adjectives classical and ballistic 
refer to the regime Af <C W <C I, where Af is the Fermi 



wavelength, W the width of the point contact, and / the 
mean free path. The width W and length L of the con- 
striction are also assumed to be much smaller than the su- 
perconducting coherence length £o = hy-p/irA® (with vf 
the Fermi velocity). It is then irrelevant for the Joseph- 
son effect whether the constriction is made out of a super- 
conductor or a normal metal. Kulik and Omel'yanchuk 
calculated the relationship between the supercurrent / 
and the phase difference 5<p of the pair potential in the 
superconducting reservoirs at opposite sides of the con- 
striction. Their zero-temperature result is 



I = ttG— sin(<W2), \6<p\ < n, 
e 



(1.1) 



where G is the normal-state conductance of the point 
contact. For a three-dimensional (3D) point contact of 
cross-sectional area S one has G = (2e 2 /h)kpS/4:ir, with 
&F the Fermi wavevector. (A 2D point contact of width 
W has G = (2e 2 / 'h)kpW / V .) These conductances, which 
are independent of the mean free path for W <C I, are 
contact conductances. The current-phase relationship 
(|1.1|) is plotted in Fig. 1 (solid curve). 

As always, the Josephson current is an odd function of 
8(j>, and periodic with period 2ir. The critical current I c = 
max 7(50) is reached at 5<fi — ir, and equals I c = nGAo/e. 
At the critical current, the I{S<j>) curve is discontinuous 
for T = 0. This discontinuity is smoothed out at finite 
temperatures. 

The supercurrent in the case of a ballistic point contact 
has a markedly different dependence on the phase than 
in the case of a tunnel junction, for which the Josephson 
effect was first proposed 01 • The T = current-phase 
relationship for a tunnel junction is given by dm 

/ = ±ttG— sm(6<b), (1.2) 
z e 



plotted as a dashed curve in Fig. 1 . There is now no dis- 
continuity in I (Sip). The critical current I c — nGAo/2e 
is reached at 5<p — n/2. In both cases, I c is proportional 
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FIG. 1 Current-phase relationship at T — of a classical bal- 
listic point contact, according to Eq. 11. H (solid curve), and 
of a tunnel junction, according to Eq. 11.21 (dashed curve). 



to G, but with different proportionality coefficients. The 
ratio I c /G is therefore not universal. 

Point contacts and tunnel junctions have in common 
that the length of the junction or weak link is much 
smaller than £o- At low temperatures, the Josephson 
effect can occur also in junctions much longer than £o 
[l3L IT3|. Consider an SNS junction (S=superconductor, 
N=normal metal) having a length Ln of the normal- 
metal region (i.e. the separation of the two SN interfaces) 
which greatly exceeds the coherence length £0 of the su- 
perconductors (while still Ln <C I)- The current-phase 
relation ship for the SNS junction at T = is a sawtooth 

HEEI 



aG 



hvp 5cj> 
eL^s tt 



\5<f>\ < TT, 



(1.3) 



with critical current I c = aGhvp/eL^ (reached at 5<p — 
tt). The numerical coefficient a is of order unity and de- 
pends on the dimensionality of the junction. (There is 
a disagreement among Refs. and 01 on the precise 
value of a.) The I c /G ratio is now a function of Ln- 
This is in contrast to the point contact and tunnel junc- 
tion, where I c /G depends only on the gap Ao in the bulk 
reservoir — not on the geometry of the junction. There is 
also a different temperature dependence. The supercur- 
rents and l|1.2J) decay when T approaches the crit- 
ical temperature T c ~ Ao/fce, wheras Eq. I|1.3|) requires 
temperatures smaller than T c by a factor £o/£n *C 1. 

The problem addressed in Ref. p| was to find the su- 
perconducting analogue of the conductance quantization 
in the normal state, which is a geometry-independent ef- 
fect. The system in which we looked for such an effect 
was a short ballistic point contact, much shorter than £ > 
and of a width comparable to Af- One would hope that 
a quantized G will lead to a "quantization" (or, more 



correctly, a discretization) of I c in units which will still 
depend on Ao, but which are independent of the prop- 
erties of the junction. The long SNS junction is not a 
likely candidate because of its J^N-dependent critical cur- 
rent. Still, if a quantum point contact is inserted in a 
long SNS junction (Ln ^ £o)j one might expect to find 
geometry- dependent quantum size effects on the Joseph- 
son current. Theoretical work on that structure was done 
by Furusaki, Takayanagi, and Tsukada 0|. Their result 
is that I c increases stepwise with increasing width for 
certain shapes of the constriction, but not generically. 
When steps do occur, the step height depends sensitively 
on the parameters of the junction. As we will show, a 
short junction (short compared to £0) is essential for a 
generic, junction-independent behavior. 



In Sec. II we describe the method that we use to cal- 
culate the Josephson current through a quantum point 
contact, which is different from the method used in Ref. 

for a classical point contact. The central result of 
that section is a relation between the supercurrent and 
the quasiparticle excitation spectrum of the Josephson 
junction. In Sec. Ill the actual calculation is presented, 
following Ref. [f|. We conclude in Sec. IV with a brief 
discussion of a possible experimental realization of the 
superconducting quantum point contact. 



The present theory is based on two approximations: 
Firstly, we assume adiabatic transport, i.e. absence of 
scattering between the transverse modes (or ID sub- 
bands) in the quantum point contact. The assumption 
of adiabaticity brings out the essential features of the ef- 
fect in the simplest and most direct way, just as it does 
in the normal state [l^. As in the normal state pcj . 
we expect the discretization of the Josephson current to 
be robust to deviations from adiabaticity which do not 
cause backscattering. In particular, although we assume 
both I ^> W and I » £ in the present analysis, we be- 
lieve that the condition I 3> W on the mean free path is 
sufficient, since scattering in the wide regions outside the 
constriction is unlikely to cause backscattering into the 
constriction. 



Secondly, in order to further simplify the problem, we 
treat the propagation of the modes in the WKB approx- 
imation. The WKB approximation breaks down if the 
Fermi energy Ep lies within Ao from the cutoff energy of 
a transverse mode (being the ID subband bottom). We 
can therefore only describe the plateau region of the dis- 
cretized Josephson current, not the transition from one 
multiple of eAo/?i to the next. Since the transition re- 
gion is smaller than the plateau region by a factor of 
order iVAo / E-p (where N is the number of occupied sub- 
bands in the constriction) , it is relatively unimportant if 
A « E F and N ~ 1. 
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II. JOSEPHSON CURRENT FROM EXCITATION 
SPECTRUM 

The theory of Kulik and Omel'yanchuk j^] on the 
Josephson effect in a classical point contact is based 
on a classical Boltzmann-type transport equation for 
the Green's functions (the Eilenberger equation |7j), 
which is not applicable to a quantum point contact. 
The present analysis is based on the fully quantum- 
mechanical Bogoliubov-de Gennes (BdG) equations for 
quasiparticle wavefunctions, into which the Green's func- 
tions can be expanded |2l| . In the present section we 
describe how the Josephson current can be obtained 
directly from the quasi-particle excitation spectrum — 
without having to calculate the Green's functions. This 
method is particularly well suited for the point contact 
Josephson junction, which has an excitation spectrum of 
a very simple form. 

The BdG equations consist of two Schrodinger equa- 
tions for electron and hole wavefunctions u(r) and v(r), 
coupled by the pair potential A(r): 



equation |21 



H 
A* 



A 

-H* 



(2.1) 



where TC = (p + eA) 2 /2m + V — Ep is the single-electron 
Hamiltonian in the field of a vector potential A(r) and 
electrostatic potential V(r). The excitation energy e is 
measured relative to the Fermi energy. Since the ma- 
trix operator in Eq. H2.ll) is hermitian, the eigenfunctions 
^ = (u, v) form a complete orthonormal set. One readily 
verifies that if (u, v) is an eigenfunction with eigenvalue 
e, then (— v*, u*) is also an eigenfunction, with eigenvalue 
— e. The complete set of eigenvalues thus lies symmetri- 
cally around zero. The excitation spectrum consists of 
all positive e. 

In a uniform system with A(r) = AqC 1 *, A(r) = 0, 
V(r) = 0, the solution of the BdG equations is 



((h 2 k 2 /2m~E F ) 2 + A 2 ) 1/2 



(2.2) 



u(r) 
v(r) 



(2e) -l/2 e -i0/2 ( £ _ h 2 k 2 /2m + Ep j 1/2 eik(2 _ 4) 



The excitation spectrum is continuous, with excitation 
gap Aq. The eigenfunctions (u, v) are plane waves char- 
acterized by a wavevector k. The coefficients of the plane 
waves are the two coherence factors of the BCS theory. 
As we will discuss in the next section, the excitation 
spectrum acquires a discrete part in the presence of a 
Josephson junction. The discrete spectrum corresponds 
to bound states in the gap (e < Ao), localized within 
£o from the junction. The Josephson current turns out 
to be essentially determined by the discrete part of the 
excitation spectrum. 

Close to the junction the pair potential is not uniform. 
To determine A(r) one has to solve the self-consistency 



A(r) 



S (r)$>*(r)u(r)[l-2/(e)], (2.5) 

£>0 



where the sum is over all states with positive eigenvalue,* 
and /(e) = [1 + exp(e/kQT)]~ 1 is the Fermi function. 
The coefficient g is minus the interaction constant of the 
BCS theory of superconductivity. At an SN interface, 
g changes abruptly (over atomic distances) from a posi- 
tive constant to zero. In contrast, the Cooper pair den- 
sity A/g goes to zero only over macroscopic distances 
(on the order of £o for a planar SN interface). Because 
of Eq. (|2.5|l . the determination of the excitation energy 
spectrum is a non-linear problem, which is further com- 
plicated by the fact that the vector potential has also 
to be determined self-consistently from the current den- 
sity, via Maxwell's equations. (The electrostatic poten- 
tial can usually be assumed to be the same as in the 
normal state.) 

Once the eigenvalue problem (|2.1|l is solved self- 
consistently, one can calculate the equilibrium average 
of a single-electron operator V by means of the formula 
(cf. Ref. |13) 

= 2 E l dY + [1 - f(e)]vVv*) . (2.6) 

Notice the reverse order, u*Vu versus vVv*, and the 
different thermal weight factors, /(e) for electrons and 
/(— e) = 1 — /(e) for holes. The prefactor of 2 accounts 
for both spin directions. (It is assumed that V does not 
couple to the spin.) We are interested in particular in 
the equilibrium current density j(r), which is given by 



j = -2- VRe(/(e)u*(p + eA)u 

oo 

+ [l-/(e)]v(p + eA)v*), 



(2.7) 



in accordance with Eq. I|2.6[) . The Josephson current / is 
the integral of j • n over the cross-section of the junction 
(with n a unit vector perpendicular to the cross-section) . 

There is an alternative (and often more convenient) 
way to arrive at the total equilibrium current / flowing 
through the Josephson junction, which is to use the fun- 
damental relation ^3 



/ = 



2e dF 



(2.8) 



between the Josephson current and the derivative of the 
free energy F with respect to the phase difference. (The 
derivative is to be taken without varying the vector po- 
tential.) To apply this relation we need to know how 



A cutoff at Hwy), with ojd the Debije frequency, has to be intro- 
duced as usual in the BCS theory. 
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to obtain F from the BdG equations. This is somewhat 
tricky, since (because of the pair potential) one can not 
express F exclusively in terms of the excitation ener- 
gies (as one can for non- interacting particles). The re- 
quired formula was derived by Bardeen et al. |23j from 
the Green's function expression for F. Here we present 
an alternative derivation, directly from the BdG equa- 
tions. 

Following De Gennes |2l|, we write 

F = U -TS, (2.9) 
U=(H) + U int , (2.10) 

Oint = - Jdr\A\ 2 /g, (2.11) 
S = -2k B £)[/ln/ + (1 - /) ln(l - /)]. (2.12) 

The energy U is the sum of the single-particle energy 
(H) [defined according to Eq. I|2.6|l ] and the interaction- 
potential energy Ui n t (which is negative, since the inter- 
action is attractive). The entropy S of the independent 
fermionic excitations can be rewritten as 



S = -2fc B ^[ln/ + (l-/)e/fc B T]. 



(2.13) 



e>0 



Using Eqs. (|2.1|l and Ij2.5|) we obtain for the single- 
particle energy the expression 



(H) = ]T[e/- e (l-/)]-2[/ int 
o 

£e /dr(|u| 2 -|v| 2 ). (2.14) 



e>0 



e>0 



The term containing the u's and v's is in fact independent 
of these eigenf unctions, as one sees from the following 
sequence of equalities: 



E 

e>0 



e f dr{\u\ 2 - |v| 2 ) = J2 I dr ( u * Hu + vHv *) 
J dr (u*Hu + v*H*v) 

(2.15) 



i T / TL 
2 I TL 



TrTL. 



Here we have made use of the completeness of the set of 
eigenf unctions (u, v) when both positive and negative e 
are included. Collecting results, the expression for the 
free energy becomes 



F = -2fc B T^ln[2cosh(e/2fc B T)] 
+ Jdr\A\ 2 /g + TrTL. 



(2.16) 



This is the electronic contribution to the free energy, 
without the energy stored in the magnetic field. Eq. 



A e 



l<p 2 




FIG. 2 Schematic drawing of a superconducting constriction 
of slowly varying width. (Taken from Ref. 



(|2.16l) agrees with Ref. [23] — except for the term TyTL, 
which is absent in their expression for the free energy.^ 

The first term in Eq. (|2.16|l (the sum over e) can be 
formally interpreted as the free energy of non-interacting 
electrons, all of one single spin, in a "semiconductor" 
with Fermi level halfway between the "conduction band" 
(positive e) and the "valence band" (negative e). This 
familiar |22l | semiconductor model of a superconductor 
appeals to intuition, but does not give the free energy 
correctly. The second term (— J7i n t) in Eq. (|2.16l) cor- 
rects for a double-counting of the interaction energy in 
the semiconductor model. The third term (Tr7i) can- 
cels a divergence at large e of the series in the first term. 
Substituting F into Eq. (|2.8|l . we obtain the required ex- 
pression for the Josephson current: 

, = -^-w*^ 

p 

- ^2k B T f de In [2 cosh(e/2fc B T)] -^7 
h J A dd(j) 



(2.17) 



where we have rewritten X) e >o as a sum over the discrete 
positive eigenvalues e p (p = 1,2,.. .), and an integration 
over the continuous spectrum with density of states p(e). 
The term Tr TL in Eq. I|2.16|l does not depend on 8(f), and 
therefore does not contribute to /. The term — U- mt does 
contribute in general. For the case of a point contact 
Josephson junction considered in this paper, however, we 
will show that this contribution (the third term in Eq. 
H2.17(l can be neglected relative to the contribution from 



t Bardeen et al. 1231 use the free energy (including the magnetic 
field contribution) to determine the self-consistent pair and vec- 
tor potentials variationally: The self-consistent A and A mini- 
mize the total free energy of the system. Since Tr H depends on 
A, it is not immediately obvious to us that it is justified to dis- 
regard this term in the variational calculation of A, as was done 
in Ref. 23] . A possible argument is that energies near the Fermi 
energy do not contribute to TrTi, because eigenvalues below i?p 
cancel those above Ep. The cancellation is not exact, however. 
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the semiconductor model. A calculation of the Josephson 
current from Eq. (|2.17[) then requires only knowledge of 
the eigenvalues. This is in contrast to a calculation based 
on Eq. (|2.7[) . which requires also the cigenfunctions. 



III. JOSEPHSON CURRENT THROUGH A QUANTUM 
POINT CONTACT 

Let us now consider an impurity-free superconduct- 
ing constriction (Fig. 2), whose cross-sectional area S(x) 
widens from S m - ln to S max S min over a length L ^> Xp 
and <C £o- The x-axis is along the constriction, and the 
area .S^ax is reached at x = ±L. We are interested in the 
case that only a small number N ~ S m i n /Xp of trans- 
verse modes can propagate through the constriction at 
energy E-p, and we assume that the propagation is adia- 
batic, i.e. without scattering between the modes. In the 
superconducting reservoirs to the left and right of the 
constriction the pair potential (of absolute value A ) has 
phase 4>i and </>2, respectively. We wish to calculate the 
current I(5(j>) which flows in equilibrium through the con- 
striction, for a given (time-independent) phase difference 
Scj) = 4>i — 4>2 S (— 7r, 7r) between the two reservoirs. The 
characterization of the reservoirs by a constant phase is 
not strictly correct. The phase of the pair potential has 
a gradient in the bulk if a current flows. The gradient is 
1 /£o when the current density equals the critical current 
density in the bulk. In our case the critical current is 
limited by the point contact, so that the gradient of the 
phase in the bulk is much smaller than 1 /£q (by a factor 
Smin/ S(x) <C 1). Since the excitation spectrum is deter- 
mined by the region within £o from the junction, one can 
safely neglect this gradient in calculating I (8(f)) from Eq. 

As we enter the constriction from, say, the left reser- 
voir, the pair potential A(r) begins to vary from the 
bulk value Aoexp(i0i), for two reasons: (a) because 
the small number of transverse modes leads to a non- 
uniform density; and (6) because of contributions from 
the N modes which have reached r from the right reser- 
voir. Non- uniformities in A due to (a) and (6) are of 
order l/N(x) ~ A F /5(x) and N/N(x) ~ S min /S(x), re- 
spectively (N(x) being the number of propagating modes 
at x). For \x\ > L we have both S(x) 3> A F and 
S(x) ^> S m m- One may therefore neglect these non- 
uniformities for |x| > L, and put 

AW_ \Aoe* ifx>Z. 

Note the difference with planar SNS junctions, where A 
recovers its bulk value only at a distance £o from the 
SN interface. The much shorter decay length for non- 
uniformities due to a constriction is a geometrical "dilu- 
tion" effect, well known in the theory of weak links 
No specific assumptions are made on the variation of A 
in the narrow part of the constriction. In particular, our 



analysis also applies to a non-superconducting constric- 
tion (A(r) = at S m i n ). Eq. i|3.1[) remains valid up to 
terms of order S m i n /S max <C 1. 

We describe the propagation of quasiparticles through 
the constriction by means of the BdG equations (|2.1|1 . 
The electrostatic potential V(r) is now the confining po- 
tential which defines the constriction. We neglect the vec- 
tor potential induced by the Josephson current. This is 
allowed if the total flux penetrating the junction is small 
compared to the flux quantum h/2e, which is generally 
the case if the constriction is small compared to £o- For 
A(r) = 0, the transverse modes \n) = 4> n (r) are eigen- 
functions of (p 2 +p 2 .)/2m + V(r), with eigenvalues E n (x). 
We expand ^(r) = Xm( u n> Y n)<t>n m t° transverse modes 
and neglect off-diagonal matrix elements (n\H\n!) and 
(n|A|n') ((. . .) denotes integration over y and z). This 
is the adiabatic approximation. The functions u„ (x) and 
v n (x) then satisfy the one-dimensional BdG equations 

(pl/2m - U„)un + A„v„ = eu„, 
-(pl/2m - U n )v n + A>„ = ev„, (3.2) 

where U n (x) = E F — E n (x) - (l/2m)(n\p 2 x \n) ss E F — 
E n (x) is the kinetic energy of motion along the constric- 
tion in mode n, and A„(x) = (n|A|n) is the projection 
of A(r) onto the n-th mode. We will consider one mode 
n < N at a time, and omit the subscript n for notational 
simplicity in most of the equations. 

In order to simplify the solution of the ID BdG equa- 
tions (|3.2H . we adopt the WKB method of Bardeen et al. 
|23|. which consists in substituting 




into Eq. I|3.2|l and neglecting second order derivatives (or 
products of first order derivatives) . One thus generalizes 
the familiar WKB method for the Schrodinger equation 
to the BdG equations, by having not only a spatially de- 
pendent wavevector k{x) — but also spatially dependent 
coherence factors exp[±i7j(x)/2]. The resulting equations 
for 7/(x) and k(x) are 

- (fr 2 /2m)krj' + e = | A| cos(r ? - </>), (3.4) 
(h 2 /2m)(k 2 -ifc') - U = i|A[ sin(?7 — (3.5) 

where A(x) = ^(x)^ 1 ^^'. In general, both 77 and k 
are complex. The WKB approximation requires that U 
changes slowly on the scale of Af , so that reflections due 
to abrupt variations in the confining potential can be ne- 
glected. Reflections (accompanied by a change in sign 
of Re k) due to spatial variations in the pair potential 



+ More precisely, the constriction should be small compared to the 
penetration depth of the magnetic field in the Josephson junc- 
tion, which length is much larger than £0 in a type II supercon- 
ductor 0. 
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are negligible provided that |A| is much smaller than 
the kinetic energy U of motion along the constriction. 
Since U > Ep — En(0), the WKB method can not treat 
the threshold regime that E-p lies within Ao from the 
cutoff energy En(Q) of the highest mode N at the nar- 
rowest point of the constriction (x = 0). The energy 
separation SE = E N+ i(0) — E N (0) ~ E F /N is much 
larger than Ao for small N, so that the threshold regime 
\Ep — -Eat (0)| < Ao consists only of small intervals in 
Fermi energy (smaller than the non-threshold intervals 
by a factor A /SE <C 1). 

For |x| > L, where A is independent of x, one has a 
constant r\ which can take on the two values rf and n h , 



C.ll 



0. 



T c,h 



arccos(e/Ao), 



(3.6) 



where <j c = 1, a h = —1. We have (j> = cf>i for x < —L and 
4> = 4>2 for x > L. The function arccost is defined such 
that arccosi 6 (0, 7r/2) for < t < 1; for t > 1, one has 
iarccosi = ln[< + (t 2 - 1) 1/2 ]. The (unnormalized) WKB 
wavefunctions (x) for \x\ > L describe an electron- 
like (e) or hole-like (h) quasiparticle with positive (+) or 
negative (— ) wavevector, 




k c ' h = (2 m /fi 2 ) 1/2 [C/ + ^'(£ 2 -Ag) 1/2 ] 1/2 - (3-8) 

The square roots are to be taken such that Re fc° :h > 0, 
Im fc c > 0, Im fc h < 0. The wavef unction IJ3. 7J) is a solu- 
tion for \x\ > L of the BdG equations up to second order 
derivatives. One can verify that for zero confining poten- 
tial the solution 13.7(1 . IJ3.8J) is equivalent to the energy 
spectrum and coherence factors given in Eqs. H2.2JI - H2.4JI . 

For e > Ao, the wavevectors fc e,h are real at |x| > 
.L, and hence the wavefunctions (J3.7JI remain properly 
bounded as \x\ — > oo. This is the continuous spectrum. 
The density of states p e ' h {e) of the electron and hole 
branches of the excitation spectrum is 



c5fc e ' h 




2tt de 



(2m/h 2 ) 1 l 2 [U + (7 c - h (e 2 - \A\ 2 ) 1 ' 2 ] 1 ' 2 - k c h . 

(3.10) 



The length L§ of the superconducting regions at opposite 
sides of the constriction is assumed to be much larger 
than £o- The wavevector change <5fc C:h is of order l/£n 
for e > A and decays as e~ 3 / 2 for e » A . It follows 
that the integral over the continuous spectrum in the 
expression (I2.17|) for the Josephson current is of order 
N(eAo/h)L/t;o, which vanishes in the limit L/£o — > 0. 

For < e < Ao, the wavevectors k e,h have a non-zero 
imaginary part at |x| > L, so that Eq. IJ3.7|) is not an 
admissible wavef unction for both x > L and x < — L. 



The two bounded WKB wavefunctions for < e < Ao 
are given at \x\ > L by 



A+V\ if x < -L, 
if x > L, 

A-^ c _ i£x<-L, 
B-V h if x > L. 



(3.11) 
(3.12) 



The transition from k h to k c on passing through the con- 
striction (associated with a change in sign of Im fc) is 
analogous to Andreev reflection at an SN interface |24j . 
Andreev reflections are to be distinguished from ordi- 
nary reflections involving a change in sign of Re fc. Or- 
dinary reflections due to the pair potential are neglected 
in the WKB approximation. By matching the wavefunc- 
tions ^± to the region |x| < L we obtain a boundary- 
value problem with a discrete energy spectrum. Since 
e < Ao <C U, we may approximate fc w ±(2mt/ /Ti 2 ) 1 / 2 
in Eq. (J3.4|) (the upper sign refers to ^>+, the lower sign 
to ^_). The boundary- value problem then becomes 

±[h 2 U{x)/2m] 1/2 7 1 '{x) + \A(x)\ cos[?7(x) - cf>(x)} = e, 

(3.13) 

7]{-L) = 0i T arccos(e/A ), 

n(+L) = 02 ± arccos(e/A ). (3.14) 

Noting that r\ is real, one deduces from Eq. IJ3.13J1 the 
inequality |r/| < (e + | A\)(h 2 U/2m)-^ 2 . Since |A| < A 
and U > E-p — E^(0), we have the limiting behavior 
\ V (L) - V (-L)\ < (£/f )(l - E N (0)/E F )-V 2 -» in the 
limit L/^q — > 0. Hence, to order L/£ the bound-state 
energy e is determined by 



arccos(e/Ao) = ±i< 



(3.15) 



independent of the precise behavior of A(r) in the con- 
striction. Since arccos(e/Ao) > 0, there is a single bound 
state per mode, with energy independent of the mode in- 
dex n < N. This iV-fold degenerate bound state at en- 
ergy A n cos((50/2) differs qualitatively from the Andreev 
levels [24j in an SNS junction with Ln 3> Co, which are 
sensitive to the mode index, to the length Ln of the junc- 
tion, and also to the pair-potential profile at the SN in- 
terfaces. This sensitivity is at the origin of the geometry- 
dependent critical current of the SNS junction which we 
discussed in the Introduction. 

We are now ready to calculate the Josephson current 
from Eq. (J2.17|) . (The alternative calculation based on 
Eq. IJ2.7JI was given in Ref. Q, and leads to the same 
final result.) The integral over the continuous part of 
the excitation spectrum does not contribute in the limit 
L/£o — » (see above). The spatial integral of |A| 2 /g 
contributes only over the region within the constriction 
(|A| being independent of S(f> in the reservoirs). Since 
A| 2 /g ~ Aq/AfCo, one can estimate 



2e 

7T 



/ 



dr\A\ 2 /g~ N(eA /h) 



L 

To' 



(3.16) 
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6<p 



FIG. 3 Current-phase relationship of a superconducting 
quantum point contact, according to Eq. I3.17H . for three dif- 
ferent temperatures (in units of Ao(0)//db): T — (solid 
curve), T = 0.1 (dashed curve), and T = 0.2 (dotted curve). 

which vanishes in the limit L/£q — *■ 0. What remains is 
the sum over the discrete spectrum in Eq. (|2.17|) . Substi- 
tuting e p = A cos(5(j)/2), p = 1, 2, . . .N, we obtain the 
final result for the Josephson current through a quantum 
point contact: 

I(Scj>) = 7v|A (T)sin(^/2)tanh f Mpcos(ty/2)) . 

(3.17) 

Since N is an integer, Eq. I|3.17[) tells us that / for 
a given value of 8<p increases stepwise as a function 
of the width of the constriction. At T = we have 
I(S4>) = N(eAo/h)sm(6<p/2), with a critical current 
I c = NeAo/h (reached at S<fi — tt). Near the critical tem- 
perature T c we have I{&4>) = N(eAl/4hk B T c ) sin(<$0), 
and the critical current is reached at 8(f) = tt/2. The 
ratio I/G (with G = 2Ne 2 /h the normal-state conduc- 
tance of the quantum point contact) does not contain N 
and is formally identical to the result for a classical point 
contact @- § 

In Fig. 3 we have plotted I(6(f>) from Eq. I|3.17|) for 
three different temperatures: T — 0,0.1,0.2 in units of 
Ao(0)/fcB = 1.76 T c . For these relatively low tempera- 
tures the T-dependence of Ao may be neglected. The 
T-dependence of the Josephson current is easily under- 
stood within the semiconductor model (cf. Sec. II), as a 
cancellation of the current due to the states in the en- 



8 The identity of the classical and quantum results for I/G should 
not be mistaken for a universal truth: It only holds within the 
approximations made in our analysis, and breaks down in par- 
ticular in the transition region between two subsequent plateaux 
of discretized supercurrent. 




FIG. 4 Top view of a possible experimental realization of a 
superconducting quantum point contact. 



ergy gap at — e p and at +e p . States above the Fermi 
energy are unoccupied at T = 0, but become occupied 
when 4fceT <^ £p- The supercurrent decays most rapidly 
near \6<fi\ ~ n, since e p is smallest there. 



IV. EXPERIMENTAL REALIZATION 

A superconducting quantum point contact (SQPC) of 
fixed width might be constructed relatively easily as a 
constriction or microbridge in a superconductor of uni- 
form composition. A test of the present theory, however, 
requires the fabrication of nanostructures with variable 
width or density. In the case of normal-state quantum 
point contacts this has been realized by depleting the 
high-mobility two-dimensional electron gas (2DEG) at 
the interface of a GaAs-AlGaAs heterostructure by ap- 
plying a negative vo ltag e to a split gate on top of the 
heterostructure 0, 0, l25| . This is an ideal system for the 
study of ballistic transport of normal electrons, because 
of the long mean free path (/ ~ 10 ^m) and Fermi wave 
length (Ap — 50 nm). 

An SQPC might be realized in the same system, if 
superconducting contacts to the 2DEG can be made [26J . 
An observation of the discretized critical current requires 
that the superconducting regions extend well into the 
constriction on either side, as illustrated schematically 
in Fig. 4 (black regions). 

A remaining problem is the mismatch in Fermi en- 
ergy that will in general exist at the interface between 
the 2DEG (where E-p ~ lOmeV) and the superconduc- 
tors (where E-p ~ leV). This mismatch, if it is abrupt, 
induces normal reflections (rather than Andreev reflec- 
tions) of quasiparticles incident on the reservoirs. The 
fabrication of sufficiently clean superconducting contacts 
to the 2DEG in a GaAs-AlGaAs heterostructure may be 
difficult, because contacting requires an alloying process. 
This problem may be circumvented by using the surface 
2DEG present as a natural inversion layer on p-InAs [27| • 
Superconducting contacts, with a shape as in Fig. 4, may 
then be evaporated directly onto the surface. Once the 
contacts have been made, a constriction of variable width 
may again be realized by means of split gates (insulated 
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from the surface 2DEG). Because of the non-linear de- 
pendence of the supercurrent on the gate voltage, such 
an SQPC has possible device applications |28j . 

Acknowledgments 

This research was supported in part by the National 
Science Foundation under Grant No. PHY89-04035. 



References 

[1] B. J. van Wees, H. van Houten, C. W. J. Beenakker, J. 

G. Williamson, L. P. Kouwenhoven, D. van der Marel, 
and C T. Foxon, Phys. Rev. Lett. 60, 848 (1988). 

[2] D. A. Wharam, T. J. Thornton, R. Newbury, M. Pepper, 

H. Ahmed, J. E. F. Frost, D. G. Hasko, D. C. Peacock, 
D. A. Ritchie, and G. A. C. Jones, J. Phys. C 21, L209 
(1988). 

[3] H. van Houten and C. W. J. Beenakker, in: Analo- 
gies in Optics and Micro-Electronics, edited by W. van 
Haeringen and D. Lenstra (Kluwer, Dordrecht, 1990); 
|cond-mat/05126Hl 

[4] B. A. Montie, E. C. Cosman, G. W. 't Hooft, M. B. van 
der Mark, and C. W. J. Beenakker, Nature 350, 594 
(1991). 

[5] C. W. J. Beenakker and H. van Houten, Phys. Rev. Lett. 
66, 3056 (1991). 

[6] H. van Houten and C. W. J. Beenakker, in: Proc. Int. 
Symp. on Analogies in Optics and Micro-Electronics, 
edited by D. Lenstra and W. van Haeringen (North- 
Holland, Amsterdam, to be published); Physica B&C 
(to be published). 

[7] For a review of the Josephson effect in weak links, see: 
K. K. Likharev, Rev. Mod. Phys. 51, 101 (1979). 

[8] I. O. Kulik and A. N. Omel'yanchuk, Fiz. Nisk. Temp. 
3, 945 (1977); 4, 296 (1978) [Sov. J. Low Temp. Phys. 3, 
459 (1977); 4, 142 (1978)]. 

[9] Yu. V. Sharvin, Zh. Eksp. Teor. Fiz. 48, 984 (1965) [Sov. 
Phys. JETP 21, 655 (1965)]. 



[10] B. D. Josephson, Phys. Lett. 1, 251 (1962). 

[11] P. W. Anderson, in: Ravello Lectures on the Many-Body 

Problem, edited by E. R. Gianello (Academic, New York, 

1963). 

[12] V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 

486 (1963); 11, 104(E) (1963). 
[13] L. G. Aslamazov, A. I. Larkin, and Yu. N. Ovchinnikov, 

Zh. Eksp. Teor. Fiz. 55, 323 (1968) [Sov. Phys. JETP 

28, 171 (1969)]. 
[14] I. O. Kulik, Zh. Eksp. Teor. Fiz. 57, 1745 (1969) [Sov. 

Phys. JETP 30, 944 (1970)]. 
[15] C. Ishii, Prog. Theor. Phys. 44, 1525 (1970). 
[16] A. V. Svidzinsky, T. N. Antsygina, and E.N. Bratus', Zh. 

Eksp. Teor. Fiz. 61, 1612 (1971) [Sov. Phys. JETP 34, 

860 (1972)]; J. Low Temp. Phys. 10, 131 (1973). 
[17] J. Bardeen and J. L. Johnson, Phys. Rev. B 5, 72 (1972). 
[18] A. Furusaki, H. Takayanagi, and M. Tsukada (preprint). 
[19] L. I. Glazman, G. B. Lesovik, D. E. KhmePnitskii, and R. 

I. Shekhter, Pis'ma Zh. Eksp. Teor. Fiz. 48, 218 (1988) 

[JETP Lett. 48, 238 (1988)] 
[20] A. Yacoby and Y. Imry, Phys. Rev. B 41, 5341 (1990). 
[21] P. G. de Gennes, Superconductivity of Metals and Alloys 

(Benjamin, New York, 1966). 
[22] M. Tinkham, Introduction to Superconductivity 

(McGraw-Hill, New York, 1975). 
[23] J. Bardeen, R. Kiimmel, A. E. Jacobs, and L. Tewordt, 

Phys. Rev. 187, 556 (1969). 
[24] A. F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964); 49, 

655 (1965) [Sov. Phys. JETP 19, 1228 (1964); 22, 455 

(1966)]. 

[25] For a review of the field of quantum transport in semi- 
conductor nanostructures, see: C. W. J. Beenakker 
and H. van Houten, Solid State Physics 44, 1 (1991); 
cond-mat/0412664 

[26] Z. Ivanov, T. Claeson, and T. Andersson, Jpn. J. Appl. 
Phys. Suppl. 26-3, 1617 (1987). 

[27] H. Takayanagi and T. Kawakami, Phys. Rev. Lett. 54, 
2449 (1985). 

[28] H. van Houten, Appl. Phys. Lett. 58, 1326 (1991). 



